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For the first time we calculate quantitatively the influence of inhomogeneities on the global 
expansion factor by averaging the Friedmann equation. In the framework of the relativistic second- 
order Zel'dovich-approximation scheme for irrotational dust we use observational results in form 
of the normalisation constant fixed by the COBE results and we check different power spectra, 
■ namely for adiabatic CDM, isocurvature CDM, HDM, WDM, Strings and Textures. We find that 

| the influence of the inhomogeneities on the global expansion factor is very small. So the error 

. in determining the age of the universe using the Hubble constant in the usual way is negligible. 

This does not imply that the effect is negligible for local astronomical measurements of the Hubble 
constant. Locally the determination of the redshift-distance relation can be strongly influenced by 
the peculiar velocity fields due to inhomogeneities. Our calculation does not consider such effects, but 
is contrained to comparing globally homogeneous and averaged inhomogeneous matter distributions. 
In addition we relate our work to previous treatments. 

m 
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00 ! I. INTRODUCTION 

Lower limits of the age of the universe are observationaly determined in many ways. Measurements of isotopic 
ratios of radioactive nuclei determine the ages of meteorites by 4.5 Gyrs Studies of the cooling of white dwarfs 

H lead to an age of our solar system of about 1 Gyr. Galactic ages could be determined to lie in a range of 12.6 to 
ON ' 19.6 Gyrs by measuring the abundance ratio of different isotopes of elements jj],f|]. Measurements of the luminosity 
of stars located at the turn-off point of the Hertzsprung Rusell diagram determine the ages of globular clusters to 
O , ; lie in a range of 12 to 18 Gyrs ^-0]. On the other hand upper limits of the age of the universe can be derived 
using cosmological models. Using a standard Friedmann-Lemaitre-Robertson- Walker (FLRW) model with vanishing 
cosmological constant (A = 0), the inverse of the Hubble parameter as measured today, Hq , provides an upper 
limit of the age of the universe (the index indicates values at the present time). Recent measurements of cepheid 
variables in the virgo cluster |||| lead to a Hubble parameter of about Ho = 80 km/(s-Mpc) (upper limit Hq 1 « 12.2 
Gyrs), leading to an age of 8.15 Gyrs for a flat universe with vanishing cosmological constant. This is far below the 
. £h ! observational lower limits cited above. There are several ways out of this dilemma: 

The first is believing in a lower value of the Hubble parameter. There are two reasons for that: firstly there exist 
other observational results [|l0| , secondly the redshift-distance relation can be influenced by the inhomogeneities, which 
will influence the Hubble constant. 

The second is to believe that the high value of the Hubble constant comes from the fact that we live in an underdensed 
region of the universe, whereas on average over the whole universe the expansion parameter is smaller |TlJ . 

The third is believing in some nonvanishing cosmological constant [^],OJ , where under some circumstances no upper 
limit can be derived (to > H^ 1 ), and the age of the universe can be about 30 Gyrs or even higher fll2]] . 

In this paper we want to investigate still another way. In the usual calculation of its age, the universe is assumed 
to be exactly isotropic and homogeneous. This might be a good approximation due to the high isotropy of the 
microwave-background-radiation, so that the FLRW description might be valid in some averaged sense. On the other 
hand the inhomogeneities are large, even on large scales, for example at a scale of ~ 10 Mpc the density constrast 
5 = bp I p might reach unity. In the early universe deviations from homogeneity and isotropy were small, but after the 
deviations became nonlinear, these inhomogeneities could influence the global expansion factor. This effect is called 
by us the backreactions of the inhomogeneities. As a result of these backreactions the value of the Hubble parameter 
cannot be taken in the usual way for a determination of the age of the universe. We will calculate quantitatively the 
effect of the backreactions and we will see how large the deviations from the usual age-determinations are. This paper 
is organized as follows: In Sec. II we present the basic equations and the averaged Friedmann equation in a general 
form. In Sec. Ill we use the results of the relativistic Zel'dovich type approxomation to second order (Russ et al. 
[Of) based on the tetrad formalism in cosmology (Kasai [pi) and calculate the backreactions using different power 
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spectra and the normalisation constant fixed by the COBE results. This paper was influenced by the pioneering paper 
from Bildhauer and Futamase we compare our results with theirs and others (Buchert, Ehlers |l6|,[l7]], Futamase 
Jl8],[l9]]) in Sec. IV. Sec. V is devoted to conclusions. 



II. THE FRIEDMANN EQUATION IN AN INHOMOGENEOUS UNIVERSE 



In this section, we summarize a general relativistic treatment to describe the non-linear evolution of an inhomo- 
geneous irrotationalQ universe [pLEOLM. The models we consider contain irrotational dust with energy-density p 
and four- velocity u M . We will neglect the curvature constant k and a possible cosmological constant A. Neglecting 
the fluid pressure and the vorticity is a reasonable assumption in a cosmological context. In comoving synchronous 
coordinates, the line element can be written in the form (Indices /i, v, ■ ■ ■ run from to 3 and indices • • • run from 
1 to 3) 



ds 2 = —c 2 dt 2 + gijdx t dx-' 



(2.1) 



and u M = (c, 0,0, 0). Then, Einstein's field equations read 
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u l . j +u k k u i . j + 3 R l 1 c 2 
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(2.2) 
(2.3) 
(2.4) 



where 3 R l j is the three-dimensional Ricci-tensor, 



-,9 m 



(2.5) 



is the extrinsic curvature, || denotes the covariant derivative with respect to the three metric gij, and an overdot (') 
denotes d/dt. We introduce the conformal factor a(t) by 



9ij = a(*) 2 7y 



(2.6) 



and we introduce the quantity V 1 -, describing the deviation from a homogeneous and isotropic expansion 



tA _ i ®W cj 1 ik ■ 



(2.7) 



Then we can write two of the Einstein equations in the following form, which we call the Friedmann equations 

^2 
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and 



oW = _4 ![ G 1 

a{t) 3c 2 ' 3 k 1 
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3a(t) 



V) 



(2.9) 



We introduce the following averaging procedure ]2S 



< A >= 



vi A ^ d3x 



(2.10) 



1 We do not consider here the effect of rotation, which might turn the effect of the inhomogeneities in the opposite direction, 
i.e. tending to increase the age of the universe. 
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where g = det<?y. V is the comoving volume of a compact domain T>(t) of the fluid [ |l 6|Jl 7| ] . V should be sufficiently 
large so that we can assume periodic boundary conditions. The scale factor ao(i) describes the expansion of this 
volume. Therefore the expansion rate of the universe is defined by 

a D (t) V a(t) k 

We then average the Friedmann equations, apply the commutation rule and neglect higher order terms (see AP- 
PENDIX |X|) to get the averaged Friedmann equations in the form 

and 

= < P > +\ < (V\? - V\V\ > . (2.13) 

CL£)(t) i>& 3 

These are the general equations for the evolution of the expansion factor of an inhomogeneous universe. They do not 
depend on a specific model of the universe. Eq. (2.13) has already been discovered by Buchert and Ehlers p^ , p7| (see 
Sec. IV). 



III. MODEL OF THE INHOMOGENOUS UNIVERSE 

We use the solution of the relativistic Zel'dovich approximation to second order (Russ et al. fll3|] ) based on the 
tetrad formalism (Kasai jL4)) to get 



a D {tf _ 8irG 
a D (t) 2 ~ !kF 



The function \l/(x) is related to the initial displacements of the particles, to first order it represents the potential of 



the density fluctuations —^'u = <$(x, t in ). Here we put ajj(t in ) = 1 and V(t in ) = 1. For a justification of eq. (3.1) see 
APPENDIX ||. We use the background relationships 

c 2 2 
Pbfan) = R n , 2 and ^ = j, r^H > ( 3 - 2 ) 

where Hq is the present value for the Hubble parameter 

2 _ a%(t ) _ 8ttG / \ _ 8ttG 



H, 



" al(t ) 3c 2 



(pb(*o) + Pcorr(^o)) = P*> (*o) (l + ^corr(^o)) ■ (3.3) 



We then get 



with 



a 2 D (t) _ Hi fa 3 D (t ) , a 2 D (t ) 



,(t) l + 5 COII (t a ) V a 3 D {t) a 2 D {t) 



,3^ 1 „ 2 ^W*o)) (3.4) 



25 f 

r(*0) = -^ 10_6/l o(l + ^in) 4 M P C- 2 ^ *'™*, m d 3 X (3.5) 



and iJo = 100/iq km/(s-Mpc). We integrate the Friedmann equation (3.4) to get the age of the universe 



_ 2 y/l + Wfa) /3 f 1 yirfi \ 

3 ffo v 2 7 Vl+^corr^o)/ " 

Since <5 C orr(io) has a negative sign the age of the inhomogeneous universe is less than the age of a corresponding- 
homogeneous one calculated with a given Hubble constant. Here, we want to estimate these differences quantitatively. 
The flat background allows a Fourier decomposition <5(x, ii n ) = V^ k <5ke jkx , so we get 
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f *> m *, m d 3 x = J2il\h\ 2 > (3-7) 
Jv k 

where |^k| 2 is the power spectrum of density fluctuations 

\5 k \ 2 =P(k). (3.8) 

What we need is to know the power spectrum at initial time ij n , where the fluctuations are still linear and just start 
to move into the nonlinear regime. We choose z- ln = 8 ^3|. The power spectrum evolves according to |m| 

P(k, Un) = * , 2 P(fc, <pr)T 2 (fc) , (3.9) 

where T{k) is a transfer function and the primordial power spectrum is assumed to be 

P(k, tpr) = Ak n . (3.10) 
The normalisation constant is fixed by the COBE results 

Acobe = ( ^) Sl^H^ (9zs) 2 . (3.11) 



V 5 J u u v T - 



L 7o 



The COBE rms-fluctuation is given to be Q rms = 9-3/xK p5|, if we take n = 1 (scale invariant primordial power 
spectrum), which is the most reasonable value |^(|. The temperature of the Microwave Background Radiation is 
T lo = 2.73K. The volume will be taken large enough so that we can convert the sum into an integral 



? - w J d3k (3 - 12) 

to get 

o ^ r 
S con (to) = -^6^ 10 ^ h2 ° {1 + Zin)2j4 / kT2 ( k ^ dk ■ ( 3 - 13 ) 

In the following we will assume ho — 0.8 and fitot = 1. For the normalisation constant one finds A w 4.35 • 10 5 Mpc 4 . 
In the following we will calculate the age of the universe using different transfer functions T{k) and power spectra 
P{k), where AkT(k) 2 = P(k, to). 



A. Adiabatic Cold Dark Matter Fluctuations 



The transfer function for adiabatic CDM fluctuations is given by 



TcDM.ad(fc) = ln(1 + ^' 34g) t 1 + 3 - 89 ? + ( 16 - 1( ?) 2 + ( 5 - 46( ?) 3 + (6.71<?) 4 ] 1/4 , (3.14) 
2.34g 



where 



ke 1 ' 2 



^cdm/iq m P c 



(3.15) 



Here 9 — p er /(1.68p 7 ) is a measure of the ratio of the energy density in relativistic particles (photons plus neutrinos) 
to that contained in photons. We will set 6 = 1, corresponding to three flavors of relativistic neutrinos plus the 
photons, and we will take ^cdm = 1- Varying f^cDM will change the result only slightly. 
The result is S COTI (to) = —2.50 • 10 -3 and the age of the universe becomes 

to « 0.9995 • . (3.16) 
6tf 
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B. Isocurvature Cold Dark Matter Fluctuations 



The transfer function for isocurvature CDM fluctuations reads 

(40g) 2 



DM,isoc 



(k) = 



1 



l + 215g+(16 9 ) 2 (l + 0.5<?) 



— + (5.6q) 8 ^ 



where 



^CDM^nMpc 1 



We again will take ^cdm = 1- 

The results is S COIt (to) = —5.63 • 10~ 4 and the age of the universe reads 

9 

t « 0.99989 ■ 



3H 



-5/4 



(3.17) 



(3.18) 



(3.19) 



C. Hot Dark Matter + Adiabatic CDM 

If we assume only one species of massive neutrinos, adiabatic fluctuations give |27| 

T„,ad(*) = exp [-0.16(feH/„) ~ {kR fv f/2\ \l + l.6q + (4.0g) 3/2 + (0.92g) s 



where 



k 



and R fv = 2.6(a,/i 2 1 )~ 1 Mpc . 



(3.20) 



(3.21) 



f^Mpc -1 

We will take f2„ = 0.3 and in the adiabatic CDM transfer function we set Ocdm = 0.7 . The total power spectrum is 
then given by P(k) = (o.3y/FjJe) + 0.7 v /P CDM (fc)) . 
The results is 6 COII (to) = —6.95 • 10 -4 and the age of the universe becomes 

2 



t w 0.99986 ■ 



3H Q 



(3.22) 



D. Warm Dark Matter Fluctuations 



Adiabatic fluctuations of warm dark matter give p7[ 

kR fw (kR fw ) 2 



(fc) w exp 



TCDM.ad, (k) , 



where 



and 



Tc 



DM, ad, 



,(*) = 1 + 1.7q + {A.'Sq) 3/2 + q 



^CDM^n M P C 1 



and 



= 0.2 ( 9cBM : dcc ) 4/3 (QcBuhD-'Mpc 



V 100 



(3.23) 



(3.24) 



(3.25) 



Here gcDM,dcc is the effective number of particle degrees of freedom when the CDM particles decoupled, values range 
from 60-300, we will set <?CDM,dec = 300 and f^cDM = 1 • 

The result is S COIY (to) — —2.90 • 10~ 3 corresponding to an age of the universe of 



to « 0.99942 ■ 



3Ho 



(3.26) 
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E. String and Texture Models 



The power spectrum of a cosmic string network evolving in a flat universe dominated by CDM is given by [E8 29 



Ak 

pft] = _ (3 27) 

[1 + a 2 k + (a 3 fc) 2 + [a A kf] [1 + (a 5 fc) 2 ] 2 ' 

where a 2 = 7.57h~ 2 , a 3 = 5.89/Jq 2 , a 4 = 1.93/i 2 and a 5 = 0.000357/ig 2 . 

This power spectum assumes that the string network is characterized by a scaling solution and that the power is 
dominated by the coherent motions of loop strings; perturbations induced by string loops are neglected p8j . 
The power spectrum of a CDM universe with perturbations seeded by textures is given by [ p0|j3l[| 

Ak 

P ( fc ) = 7 757V > ( 3 - 28 ) 

[1 + {ak + (/%) 3 / 2 + [ikfY] 1 

with v = 1.2, a = 19. Ah^ 2 , (3 — 6.6/iq 2 , and 7 = 3.0ft.Q 2 , where we still set Qq = 1. 

Although the non-Gaussian nature of the string and texture models means that the power spectrum does not provide 
a full description of the density field even in the linear regime, the power spectrum is sill a well-defined quantity, and 
its meaningful to compare it to observations |^8| and give excellent fits to it (3^] . 

The result for the string and texture models are almost exactly the same and give 6 COTr (to) = —8.88 • 10 -3 . The age 
of the universe in this case becomes 

t « 0.9982 • -|- . (3.29) 
dun 



IV. COMPARISON WITH PREVIOUS WORKS 



A. Newtonian treatment 



For a comparison with the Newtonian treatment by Buchert and Ehlers Jl6|,[T7j we have to identify their Vv with 



our it\, and their Vu with our V %. Their result 



ap 



AttGM 1 . , ., 
— V + 3 <{U fe) 



> 



(4.1) 



is found to agree exactly with ours in eq. (2.13), except for the derivatives, which are covariant in our case. Eq. (23) in 



is an extension to describe a globaly anisotropic universe. They concluded, that their eq. (9) in |16| and eq. (19) 
11711 must also hold in general relativity, because they were derived by a vera ging the Raychaudhu ri eq uation. In 



our treatment the Raychaudhuri equation can be derived by combining eq. (2.2) and the trace of eq. (2.4), replacing 



u' d = (1/3)6^ - 
in our treatment. 



o % a. So it is possible to recover all their results, exept for the vorticity, which we assume to vanish 



B. Futamase's approximation scheme 

Futamase Jl8| , ^9[ calculated the backreactions based on his approximation scheme, where he introduced two small 
parameters representing the amplitude of the metric fluctuations and the ratio between the scale of the variation 
of this metric fluctuations and the scale of a(t) and the background metric. Then in [ff8| he used a cosmological 
post-Newtonian approximation. In | fl9[ he employed the 3 + 1 splitting of space time, then the Isaacson averaging 
J32l | is performed on the background spatial hypersurface. His results in his eq. (3.16) or in eq. (68) [[l9| are of 
the same order as ours, but the factors are different. There are several reasons for that discrepancy: firstly in his 
approximation sheme he neglected some terms, we don't use such an approximation. Secondly he did not introduce a 
scale factor ajj(t) defined by the expansion of a comoving volume. He introduced the conformal factor a(t), then he 
rescaled this scale factor by neglecting terms like < h k k >. Thirdly his averaging process in [Q is not defined using 
the square-root of the real metric under the integral, rather he used the square-root of the background metric, which 
is essentially unity for a flat background. He also used this averaging process at the end of (l9| to recover the results 
of@. 
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C. Bildhauer and Futamase 



Bildhauer and Futamase ||15(] calculated the backreactions of the inhomogeneities based on the work of Futamase 
and the Newtonian Zel'dovich approximation (Buchert p3[). Their result (eq. (25), see also eq. (84) in [Q) reads 

With b = (Jcorr(io) 

<W*o) - S!l0- 6 ^(1 + Z in ) 4 ^ < I^P > , (4.2) 

where we want to indicate a typing error (Mi defined in their eq. (25) is not the same as in their eq. (22), the factor 
k 2 1 jjL 1 is incorporated into Mi). With pT 2 < \U\ 2 >= [sT 2 < |Vsi n | 2 >=< \J" m ^, m > this is of the same order as 
our result, only the factor is different. The reasons are the same as those in the previous subsection. Another error 
was found by Futamase fL9|| : Mi = 57ir 3 should read Mi = 57/8, the mistake is coming from the wrong integration 
region [0,27r] instead of [0,d]. They derived at the conclusion that the underestimation of the age of the universe 
is approximately 30 percent, which is not correct since they just assumed (5 CO rr(io) to be of order unity instead of 
calculating it quantitatively as we did here. 



V. CONCLUDING REMARKS 



We have calculated quantitatively the influence of the inhomogeneities on the global expansion factor of a flat 
universe with vanishing cosmological constant in the framework of a Zel'dovich type relativistic approximation scheme 
using the results from COBE. The first result is that the backreactions act as an additional energy density, which 
is proportional to a~jj 2 , so we can interpret the averaged expansion as Friedmannian with a small positive spatial 
curvature. The second result is that this influence is very small. As a consequence of this the modification of 
the age of the universe calculated in the usual way (i.e. assuming a homogeneous universe) with a given Hubble 
constant is negligible. In all models considered here relative differences were less than as 2 • 1CP 3 . This does not 
imply that the inhomogeneities are negligible for local astronomical measurements of the Hubble constant. Locally 
the determination of the redshift-distance relation can be strongly influenced by the peculiar velocity fields due to 
inhomogeneities (d = Hz + 0(2)). Our calculation does not consider such effects, but is contrained to comparing 
globally homoeneous and averaged inhomogeneous matter distributions. Calculating the modification of the redshift- 
distance relation will be the subject of future investigations. As a result the age problem of the universe that arises 
in high density models can only be solved either with a lower Hubble constant, with a nonzero cosmological constant 
or with a reduced age of globular clusters. 
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APPENDIX A: COMMUTATION RULE 



The time derivative of an averaged quantity reads: 



l<A>=-^<A>+± I (A^g + AyTg)d 3 x. (AH 



V 



This leads to the commutation rule |22,lq,l'i 



where 



j t <A>-<A>=-<6><A> + <A6>, (A2) 



and <0>=3 ^W. (A3) 

y/g a D(t) 
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To convert eqs. (2.8) and ( fe.gQ to the eqs. ( 2.1 2| ) and ( [2.13] ) we used 

d 



k ^2 



— < vi > - < vi >=< (v%y >-<v%> 
iii 

and neglected the term < V \ > 2 , because it is a higher order quantity. 



APPENDIX B: MODEL OF THE INHOMOGENEOUS UNIVERSE 



The result of the relativistic Zel'dovich approximation to second order [|13[ is the following metric tensor: 

20 



1 



9c 2 t? ^)S ij +2a D (t)^ ttJ 



apjt) 
cH 2 

c t in 

■a 2 D (t) 



20 40 T T 10 k t r 

— 3 - g-**,ij + g-* , ** > fc<5« 

19 T fe T 12 T fe T 3 // fe 



where we set a = —50/81. Since their difference is only of second order, we could replace a(t) by an 
determinant we only need to first order: 



V7 = l 



1 , . T k 



3 c 2 *? 



This leads to 



< A R > 



2 2+2 \r , 
a D C l in V m JV 



1 / 400 



9 < m 



600 



an 



81 



) Id 3 a; 



and 



< (v\f - V k m V" k l >= 



9a D t?Vi 



in I'm JV 



The averaged Friedmann equations read: 



a 2 D _ 8ttG 1 



20 
27a 2 D 



4" 



27od 



f c*. m ^i 2 _ vF fe 



m ,fe 



af^t 2 , V243 



300 



\F m vI/ 



200 
i 

243 



U 3 x 



and 



an _ 47rG 
3c 2 



k ,Tr> m ^ ,j3 



(Fx 



Integration and the assumption of periodic boundary conditions lead to 



^' m m d a x = . 



The Fourier transformation <5(x, t in ) = ^ k <5ke lkx together with J v exp l ( k k ) x d 3 x = <5kk' leads to 



fd 3 x 



and 



f W™ n d 3 x = - [ ^' m ^, m <fx . (B9) 
Jv Jv 

)ws: 

< P >= -^77 y p(tin)v / 5(Urf 3 ^ = — y = -^TPb(t) , (BIO) 

where we used the conservation law p{t- ln )\/ g(t ln ) = p{t)y/g(t). Note that even in a case where the averaged second 
scalar invariant would not vanish our treatment would still be consistent, since in every case 



The averaged density is treated as follows: 

8ttG 8ttG 1 / ., — 87,-G M SttG 



— ( tlh.\ - 2^5. ( 

dt \a 2 D J a D \a D a 2 D 



■2 ;- 2 ^ r^- J) (Bn) 



is satisfied. 
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